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Abstract. We study the following class of scalar hyperbolic conservation laws with discontin- 
uous fluxes: 

dtp + d^F{x,p) = 0. (0.1) 
The main feature of such a conservation law is the discontinuity of the flux function in the 
space variable x. Kruzkov's approach for the -contraction does not apply since it requires the 
Lipschitz continuity of the flux function; and entropy solutions even for the Riemann problem 
are not unique under the classical entropy conditions. On the other hand, it is known that, 
in statistical mechanics, some microscopic interacting particle systems with discontinuous speed 
parameter \{x), in the hydrodynamic limit, formally lead to scalar hyperbolic conservation laws 
with discontinuous fluxes of the form: 

dtp + d:,(\{x)h{p)) =0. (0.2) 

The natural question arises which entropy solutions the hydrodynamic limit selects, thereby lead- 
ing to a suitable, physical relevant notion of entropy solutions of this class of conservation laws. 
This paper is a first step and provides an answer to this question for a family of discontinuous 
flux functions. In particular, we identify the entropy condition for HO. Il l and proceed to show the 
well-posedness by combining our existence result with a uniqueness result of Audusse-Perthame 
(2005) for the family of flux functions; we establish a compactness framework for the hydrody- 
namic limit of large particle systems and the convergence of other approximate solutions to 110. Il l , 
which is based on the notion and reduction of measure-valued entropy solutions; and we finally 
establish the hydrodynamic limit for a ZRP with discontinuous speed-parameter governed by an 
entropy solution to 1 10.21 1. 



1. Introduction 

We are concerned with the foUowing class of scalar hyperbolic conservation laws with discon- 
tinuous fluxes: 

dtp + d^F{x,pit,x))=0 (1.1) 

and with initial data: 

p\t=o=PQ{x), (1.2) 
where F{-,p) is continuous except on a set of measure zero. 

The main feature of (jl.ip is the discontinuity of the flux function in the space variable x. 
This feature causes new important difficulties in conservation laws. Kruzkov's approach in |18j 
for the i^-contraction does not apply; entropy solutions even for the Riemann problem of p.ip 
are not unique under the classical entropy conditions; several admissibility criteria have been 
proposed in [H [3l HI [151 ttZ] and the references cited therein. In particular, a uniqueness theorem 
was established in Baiti-Jenssen |3] when F{x, •) is monotone and Audusse-Perthame [T] for more 
general flux functions that especially include non-monotone functions F{x, •) in (jl.ip under their 
notion. However, the existence of entropy solutions for the non-monotone case under the notion 
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of Audusse-Perthame [IJ has not been established, and the entropy conditions proposed in the 
hterature in general are not equivalent. 

On the other hand, in statistical mechanics, some microscopic interacting particle systems with 
discontinuous speed parameter A(a;), in the hydrodynamic limit, formally lead to scalar hyperbolic 
conservation laws with discontinuous flux of the form 



and with initial data (|1.2p . where A(a;) is continuous except on a set of measure zero and h{p) 
is Lipschitz continuous. Equation (|1.3p is equivalent to the following 2x2 hyperbolic system of 
conservation laws: 



In particular, when h(p) is not strictly monotone, system (jl.4p is nonstrictly hyperbolic, one of the 
main difhculties in conservation laws (cf. [5l[7]). The natural question is which entropy solution 
the hydrodynamic limit selects, thereby leading to a suitable, physical relevant notion of entropy 
solutions of this class of conservation laws. This paper is a first step and provides an answer to this 
question for a family of discontinuous flux functions via an interacting particle system, namely, 
the attractive zero range process (ZRP). This ZRP leads to a conservation law of the form (|1.3p 
with A(a;) > and h{p) being monotone in p, and its hydrodynamic limit naturally gives rise to 
an entropy condition of the type described in [1] |^ in the formal level. 

Motivated by the hydrodynamic limit of the ZRP, in this paper, wc adopt the notion of entropy 
solutions for a class of conservation laws with discontinuous flux functions, including the non- 
monotone case in the sense of Audusse-Perthame , and establish the existence of such an entropy 
solution via the method of compensated compactness in Section 3. This completes the well- 
posedness by combining a uniqueness result established in [T] for this class of conservation laws 
under the notion of entropy solutions. 

In order to establish the hydrodynamic limit of large particle systems and the convergence of 
other approximate solutions to (jl.ip rigorously, we establish another compactness framework for 
(|l.ip - (|1.2p in Section 2. This mathematical framework is based on the notion and reduction of 
measure- valued entropy solutions developed in Section 2, which is also applied for another proof 
of the existence of entropy solutions for the monotone case in Section 3. 

In Section 4, we establish the hydrodynamic limit for a ZRP with discontinuous speed-parameter 
A(a;) governed by the unique entropy solution of the Cauchy problem (|1.2p - (|1.3p . 



In this section, we first develop the notion of measure-valued entropy solutions and establish 
their reduction to entropy solutions in L°° (provided that they exist) of the Cauchy problem 
(dHll-lIOl) satisfying 

(HI) F{x,p) is continuous at all points of (IiyA/") x H with J\f a closed set of measure zero; 
(H2) 3 continuous functions f,g such that, for any a: G K, and large p, f{p) < \F{x,p)\ < g{p) 

with f{p) > and /(±oo) — oo; 
(H3) There exists a function Pm{x) from R to R and a constant Mq such that, for x E M\J\f, 

F(x, p) is a locally Lipschitz, one to one function from (— cxd, p^] and [pm, oo) to [Mq, oo) 

(or (-00, Afo]) with F{x,pm{x)) = Mq; 



dtp + d^{\{x)h{p))=Q 



(1.3) 




(1.4) 



2. Notion and Reduction of measure-valued entropy solutions 



or 



(113') For X G HyA/", F{x, •) is a locally Lipschitz, one to one function from R to R. 
One example of the flux function satisfying (H1)-(H2) and (113) or (113') is 

F{x,p)^X{x)h{p), 



(2.1) 
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where X{x) is continuous in a; € H with < Ai < X{x) < A2 < oo for some constants Ai and 
A2, except on a closed set Af of measure zero, h{p) is locahy Lipschitz and is either monotone or 
convex (or concave) with h{pm) — for some pm in which case Mq = 0. 

It is easy to check that, if the flux function F{x,p) satisfies (H1)-(H3), then, for any constant 
a G [A/0,00) (or a S (— oo,Mo]), there are two steady-state solutions from II to [pm{x),oo) 
and from R to (— 00, p„i(a;)] of (jl.ip such that 

F{x,mt{x))=a. (2.2) 

In the case (H1)-(H2) and (H3'), m'^{x) — m~{x) which is even simpler. 

2.1. Notion of measure-valued entropy solutions. First, the notion of entropy solutions in 
L°° introduced in Audussc-Pcrthame 1 and Baiti-Jcnssen [3j can be further formulated into the 
following. 

Definition 2.1 (Notion of entropy solutions in L°°). We say that an L°° function p : := 
R-)_ X R i—f R is an entropy solution of p.ip - (ll.2p provided that, for each a G [Mo,oo) (or 
a e (— oOjMq]^ and the corresponding two steady-state solutions m^{x) of (jl.ip . 

\p(t, x) — m^{x) \ dtJ + sgn{p{t, x) — m^{x))(^F{x, p{t, x)) — a) d^J^ dtdx 

\pQ{x)~m^{x)\J{Q,x)dx>Q (2.3) 

for any test function J : 1-^ 1R+ . 

It is easy to see that any entropy solution is a weak solution of (|l.ip - (|1.2p by choosing a such 
that m^ix) > IIpIIl^ and m~{x) < —\\p\\l'^, respectively, for a.e. x G R. 

From the uniqueness argument in Audusse-Perthame !T (also see [^j), one can deduce that, for 
any L > 0, 



lim / \p{t,x) - po{x)\dx = 0. (2.4) 

*^°J\x\<L 

Following the notion of entropy solutions, we introduce the corresponding notion of measure- 
valued entropy solutions. We denote by ^(IR) the set of probability measures on R. 

Definition 2.2 (Notion of measure-valued entropy solutions). We say that a measurable map 

TT : ^ r{R) 

is a measure-valued entropy solution of (ll.ip - (|1.2p provided that {Tro^x',k) = po{x) for a.e. x eR 
and, for each a € [Mo,oo) (or a G (— oo,Mo]j and the corresponding two steady-state solutions 
m±(x) of 

{{'^t,x; \k - m^{x)\)dtJ + {-Kt^x] sgn{k - m^{x)) {F{x, k) - a))dxJ) dxdt 

Pf){x)-mt{x)\J{Q,x)dx>{) (2.5) 

for any test function J : R'^ 1-^ R+ . 

If a measure- valued entropy solution T:t,x{k) is a Dirac mass with the associated profile p{t,x), 
i.e. TTt xik) = (5p(t x){k), then p{t,x) is an entropy solution of p.ip - (|1.2p . which is unique as shown 

Note that, when the flux function F(x,p) is locally Lipschitz in both variables {x,p), one can 
use the Kruzkov entropy inequality, instead of (|2.5p , to formulate the following notion of measure- 
valued solutions: 

dt {TTt,x; \k - c|) -f {TTt,x; sgn(fc - c) {F{x, k) - F{x, c))) + (^t,.; sgn(fc - c)d,F{x, c)) < (2.6) 
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in the sense of distributions and to establish their reduction as in DiPerna [12j . One of the new 
features in our formulation (|2.5p in Definition 2.2 is that the constant c in (j2.6p is replaced by 
the steady-state solutions mj(x) so that the additional third term in (|2.6p vanishes, as in [Bill) 
and thereby allows the discontinuity of the flux functions on a closed set of measure zero for 
measure- valued entropy solutions. 

2.2. Reduction of measure- valued entropy solutions. In this section we first establish the 
reduction of measure-valued entropy solutions of (|l.ip - (jl.2p and prove that any measure-valued 
entropy solution TTt^x{k) in the sense of Definition 2.2 is the Dirac solution such that the associated 
profile p{t, x) is an entropy solution in the sense of Definition 2.1. That is, our goal is to establish 
that, when iro.xik) = (5p^(a.)(fc), 

7rt,a;(fc) = <5p(t,2;)(fc), (2.7) 

where p : ^ R is the unique entropy solution determined by (|2.3p . The reduction proof is 
achieved by two theorems. We start with the first theorem. 

Theorem 2.1. Assume p : ^ R is the unique entropy solution of (|l.ip - (|1.2p with initial 
data pq G L°°(R). Assume that there exists a measure-valued entropy solution n : R'^ V{R) 
of (|l.ip in the sense of Definition YI.'J] with nt.x having a fixed compact support for a.e. (t,x) and 
T^o.xik) = (5pQ(a;)(fc) for a.e. x G R. Then 

J ( (TTt,,; \k - p{t, x)\) dtJ + {TTty, sgn{k - pit, x))iFix, k) - F{x, pit, x)))) d^J) dxdt > (2.8) 

for any test function J : R^ R+ . 

Proof. The proof is divided into six steps. 
Step 1. We first rewrite 

E := dt{-Kt,x; \k - pit,x)\) dx{'nt,x]Sgnik- pit, x)) {Fix,k) - Fix, pit, x)))) 

in the entropy inequality (|2.8p . We notice the following: 

• Under the assumption (H3'), Fix, pit, x)) is continuous in x a.e. Then we can define a 
function pis,y,x) for a.e. is,y,x) G R+ x R^ such that, for fixed (s,y), 

Fix,pis,y,x)) := Fix,mF(^y^p{s,y))ix)) = F(y,p(s,y)), (2.9) 

where the last equality follows from (12. 2p . Thus, we define 

Pis,y) -.^ Fiy,pis,y)) so that pis,y,x)^mp(s,y)ix). 

• For the case (H3), we define the sign of the difference between the tilda function and Pmiu) 
to be the same as the sign of the corresponding solution. Since Pmiy) is the minimum (or 
maximum) point of the flux function with Fiy, pmiy)) = Mq, then, for 

pis,y,x) := m+(^ ,^j(a;)sgn+(p(s,y) - p„(y)) + m"^^ (x)sgn_ (p(s, y) - Pmiv)), (2.10) 

we have as in (I2.9p 

Fix,pis,y,x)) := i^(a;, to+ ^ ,^)(a:)sgn+(p(s, y) - p„iiy)) + m~^^ y~^ix)sgn_ipis,y) - p,niy)) 
= Fiy, pis, y)) = Pis,y). 
With these notations, we set 

E -.^ dt{nt^x;\k~ pis,y,x)\) + dx{TTt,x;sgnik ~ pis,y, x)){Fix, k) -/3(s,y))). (2.11) 
Then, to obtain the inequality £' < 0, it suffices to show that lima;^^ E — E. 
Step 2. We now show that 

pis, y, x) "-^^ pis, y, y) = pis, y) for a.e. (s, y) G R+. (2.12) 
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For the case (H3'), since the flux function is continuous outside a negligible set Af, then, for 

F{x, p{s, y, y)) F{y, p{s, y, y)). 
On the other hand, we have F{y, p{s, y, y)) = F{x, p{s, y, x)). Therefore, we have 

F{x, /5(s, y, x)) - F{x, p{s, y, y)) —i 0, 

and l|2.12p is a consequence of the fact that F{x^ ■) is a one to one function. The case (H3) is clear 
from the definition of p{s,y,x) in (|2.10p . 

Step 3. With Steps 1-2, to achieve inequality (|2.8p . it suffices by choosing a = (3{s,y) in (|2.5p 
to show the following inequality: 

lim / dtJ{t,x)HT{t ~ s)Hui{x - y){TTt,x]\k - p{s,y,x)\) dtdxdsdy 



+ lim / J(Q,x)Hr{—s)Huj{x — y)\po{x) — p{s,y,x)\dxdsdy > (2-13) 

for any test function J G C^{R'^) and verify that p{s, y, x) can be replaced by p{t, x) in the limit 
as r, ^ 0. Here the two families of functions Hr^ e C^(R) are defined as 

H(z) = -H{-) and HJz)^-H(-) forr,w>0, 
T r uo uj 

for a positive, compactly supported function H e C^(]R) and a positive function H e C^(]R) 

with compact support in (—2, —1) such that Jj^H(z)dz — 1 and Jj_2 _i) H{z)dz — 1. 

This can be easily seen by first choosing the test function in (|2.5p as J(t, x)Hr{t—s)H^{x—y) > 

for fixed {s,y) and then integrating with respect to (s,y). We now estimate the three terms of 

(|2.13p in Steps 4-6, respectively. 

Step 4- We show that, as r, ^ 0, the first term converges to 

dtJ{t,x) {-Kty, \k - p{t, x)\) dtdx. 
Observe that 

dtJ{t,x)Hr{t - s)H^{x - y) {-Kt^x] \k - p{s,y,x)\) dtdxdsdy 



-J dtJ{t,x)Hr{t - s)H^{x ~ y){T:t,x] \k - p{s , y , y)\) dtdxdsdy 
< / \dtJ{t,x)\Hr{t- s){ / H^{x ~ y)\p{s,y,x) - p{s,y,y)\dx) dtdsdy ^ as 0, (2.14) 



by the dominated convergence theorem and the fact that J Hi^(x — y)\p(s, y, x) — p{s, y,y) \ dx —> 



as w ^ for a.e. {s,y) e 11+ since p{s,y,x) p{s,y,y) = p{s,y) by (|2.12p . Then, to find the 



limit of the first term of (|2.13p . it suffices to compute the limit of 

dtJ{t,x)Hr{t - s)H^{x - y){'Kt,x; \k - p(s,y)|) dtdxdsdy. (2-15) 
Thus, it suffices to show that p{s,y) can be replaced by p{t,x) in ()2.15p . i.e., as r, a; ^ 0, 
dtJ{t, x)Hr{t — s)H^[x — y)\p{t, x) — p{s, y) | dtdxdsdy 

dtJ{t, x)H{-r)H{-z)\p{t, x) - p{t + rr, a; + uoz) I dtdxdrdz 0. (2.16) 
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This is guaranteed by the fact that 



hm / \p{t, x) — p{t + rr, X + Luz)\ dtdx ^ 0, 

and the dominated convergence theorem since all the functions involved are bounded. This implies 
that, in (|2.15p . we can indeed replace p{s,y) by p{t,x) to arrive at the result. 
Step 5. We show that the second term of (|2.13p converges to 



dxJ{t, x){'Kt,x]Sgn {k — p{t, x)) [F{x, k) — f3{t, .x))) dtdx as t, w ^ 0. 

The hypothesis (H2) on F{x, p) implies 

sgn (fc - p{s, y, x)) {F{x, k) - f3{s, y)) - sgn (fc - p{s, y, y)) {F{x, k) - F{x, p{s, y, y))) 

= sgn (fc - p{s, y, x)) {F{x, k) - F{x, p{s, y, x))) - sgn (k - p{s, y)) {F{x, k) - F{x, p{s, y))) 
< C\p{s,y,x) - p{s,y)\. 

Integrating the last expression with respect to x against the function H^^ (x — y) yields its conver- 
gence to zero by Step 2 as w — > 0. Therefore, the limit of the second term of (|2.13p is the same as 
the limit of 

dxJ{t,x)Hr{t ~ s)H^{x - y){'Kt,x;sgn{k-p{s,y))[F{x,k) - F{x, p{s,y)))) dtdxdsdy, 

and it suffices to prove that, as r, w — > 0, 

/ dxJ{t,x)Hr{t- s)H^{x - y){nt.x; \ sgn {k - p{s , y)) {F{x,k) ~ F{x,p{s,y))) 

—sgn {k — p{t, x j) [F{x, k) — F{x, p{t, x))) |) dtdxdsdy 0. 

Using the Lipschitz property and fact (|2.16p . we achieve the result for the second term of (|2.13p . 

Step 6. We now show that the third term of (|2.13p converges to zero if r, w — > 0. Note that 



J{0,x)Hr{-s)Huj{x - y) \po{x) - p{s,y,x)\ - \po{x) - p{s,y,y)\ 



dxdsdy 



< / J{0,x)Hr{-s)H^{x ~ y)\p{s,y,x) - p{s,y,y)\dxdsdy. 



For the same reason as in the first and the second term of (|2.13p , the right hand side converges to 
zero if r, w 0. We therefore next compute the limit as t, w ^ of 

J{0,x)Hr{-s)H^{x ~ y)\po{x) - p{s,y)\ dxdsdy. 

As before, ImiT^uj^o J J(Q,x)Hr{-~s)H^{x — y)\p(s,x) — p{s , y)\ dxdsdy = 0. Therefore, the next 
goal is to compute the limit of 



J(0, x)Hr{~s)H^{x — y)\pQ{x) — p{s, x) \ dxdsdy = J J(0, x)H{~r)\po{x) — p{Tr, x) \ dxdr. 

(2.17) 

Since all the functions are bounded and swpp H C (—2, —1), by the dominated convergence theorem 
and property (|2.4p of the unique entropy solution p(t, x), this converges to as r and thereby 
p.l7p converges to 0. 

With Steps 3-6 and by (|2.13p . we complete the proof. □ 

Then Theorem 2.1 yields the L^-contraction between the measure-valued entropy solution irt^x 
and the unique entropy solution p{t,x) of (|l.ip - l|1.2p . 
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Theorem 2.2 (L^-contraction). The function J{Trt,x', \k — p{t, x)\) dx is non-increasing int > 0, 
which implies iTt.xik) = 6p(t.x){k) when TTo^x{k) = (5po(a;)(fc) for a.e. x G R. 

Proof. In expression (|2.8[) . we choose the test function as the product test function Jj{t)H{x), with 
Jj{t) converging to the indicator function Ij^^ j^jC^) as j ^ oo for t2 > ti > 0. Then (|2.8p is equal 
to 



H{x){'Kti,x{k)] \k - p(ti,x)\) dx- J H{x){'Kt2,x{k); \k - p{t2,x)\) dx 

H'{x){'Kt,x{k)\^g^{k - p{t,x)) {F{x,k)- (3{t,x)))dxdt > 0. (2.18) 



In (|2.18p . we choose 



for X € C!^{—2, 2) with x{^) — 1 when x E [—1, 1] and x{^) ^ 0- Letting N ^ oo first and 7^0 
then yields that, for any t2 > ti > 0, 



{T^t2,x]\k - p(t2-,x)\) dx ~ j {'Ktiy,\k- p{ti,x)\)dx < 0. 
In particular, when ^2 = ^ > 0,ti ^ 0, then nQ^x{k) — '5po(x)(fc) implies 

{'n-ty, \k- p{t,x)\)dx < 

so that ITt.xik) = 6p{t.x){k) for any t > 0. □ 

3. Existence of entropy solutions 

In this section, we establish the existence of entropy solutions (|l.ip - (|1.2p in the sense of Defi- 
nition 2.1, as required for the reduction of measure-valued entropy solutions. More precisely, for 
each fixed e > 0, denotes the unique Kruzkov solution of ()l.ip - ()1.2p in the sense p.3p . where 
the flux function depends smoothly on the space variable x; then it is shown that the sequence 
converges to an entropy solution of (jl.ip - (|1.2p . 

3.1. Existence of entropy solutions when F is smooth. Define Fe{x,p) the standard molli- 
fication of i^(a;, p) in a; G R: 

F,{x,p) {F{-,p)*0^){x) F{x,p) a.e. as e ^ 0, (3.1) 

with 6'^{x) := ^(f ), 6{x) > 0, snpp9{x) C [-1, 1], and j\ 9{x)dx = 1. For fixed e > 0, consider 
the following Cauchy problem: 

/ dtp + dxF,{x,p) = 0, , . 

\ p\t=o = Po{x) > 0. ^■'■^> 

Kruzkov's result in [TH] indicates that there exists a unique solution of (|3.2p satisfying the 
Kruzkov entropy inequality: 

dt\p'{t, x)~c\+ dx{sgn{p'{t, x) - c){F,{x, p^{t, x)) - F,{x, c))) + {sgn{p'{t, x) ~ c)dxF,{x, c) < 

(3.3) 

in the sense of distributions. We now show that the entropy solution also satisfies (|2.3p . 

Proposition 3.1. Let p'^{t,x) be a solution of the Cauchy problem (|3.2p satisfying the Kruzkov 
entropy inequality p.3p . Then p'^(t,x) also satisfies the entropy inequality (|2.3p with steady-state 
solutions m 



m 
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Proof. In (|3.3p . we choose the constant c = m'^^{y) for any a £ [A/o,oo) (or a € (— oo,Aio]), 
integrate against a test function Ji^{t,x,y) :— J{t, ^^^^)Hi^{x — y) with as in the proof of 
Theorem 12. 1[ integrate by parts in the term involving H^{x — y) with respect to dy, and observe 
that {dx + dy)Juj{t,x,y) = dxJ{t, ^^)H^{x — y) to obtain from (|3.3p that 



x) - m''^'^{y) \ H^{x - y)dtJ dtdxdy 
+ / sgn(p^(t,x)-m^'±(y)) {F,{x, {t,x)) - F,{x,m'^^{y)))H^{x - y)dyJ dtdxdy 
sgn {p''{t,x) - m^'^(y)) dyF^{x,rrf^'^ {y))H^{x - y) J dtdxdy 
sgn {p''{t,x) - m^'=^(y)) dxFs{x,m%{y))H^{x - ?/)J dtdxdy 



+ J \p'iO,x)~m^^iy)\HUx-y)JiO,^^)dxdy>0, (3.4) 

where we have used that 

{F,{x, p^t,x)) - F,{x,m'^^{y)))HUx - y)J dy{sgn{p^t,x) - m'^^{y)))dxdt ^ 0. 

As in the proof of Theorem 12.11 we can replace p'^{t,x) by p'^{t,y) in the first term as w — > 
and replace p^(0, x) by p^(0, y) in the last term as w ^ 0. 
The second term is equal to 

sgn {p%t,x) - m^'=^(y)) {F,{x, p%t, x)) - F,{x,ml;'^{y)))H^{x - y)dyJ dtdxdy. 

By the hypothesis (H3) on the flux function, Fg{-,p) is a Lipschitz function from {—oo, p„i] and 
[pm.oo) to [Mq,oo) (or (— oo,Mo]), which implies 

|sgn {p'{t, x) ~ p%t, y)) {F{x, p%t, x)) - F{x, p%t, y))) \ < C\p'{t, x) - p'{t, y)\. 

One can show in a similar way as in the first term that 

lim / \p'^(t,x) — p'^{t,y)\H^{x — y)dyJ dtdxdy — 0. 

This means that, in the second term of (|3.4[) . one can replace F{x, p'^{t, x)) by F{x, p^(i, y)). Since 
Fg is also a smooth function with respect to the first variable, the second term converges to 

sgn{p'{t,y) - m^'=^(y)) {F^{y, p^{t,y)) - a)dyJ{t,y) dtdy. 

In the third and fourth term in (|3.4p . for z € R, we have 



lim / {dx+dy)F^{x,m''^^{y))ujH^{ujz)J{t,y + \ujz)dz 

cjj— >0 J A 

= lim l{d.,^dy)F,{x,Trf^^{y))H{z)J(t,y+\ujz)dz 

= J(<, y)dyF,{y, m^-±(y)) = J(t, y)dy(x = 0. 

With these results, as uj^ 0, inequality (|3.4p becomes (|2.3p for Fe{x,p) = {F{-,p) * 0'^){x) with 
steady-state solutions — rrf^^. □ 

Thus we conclude the existence of an entropy solution p^{t,x) in the sense of Definition 2.1 for 
each Fp with fixed e > 0. 
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Remark 3.1. Notice that the sequence of approximate entropy solutions converges to a measure- 
valued entropy solution as £ — s- 0: First, since po e L°°, we find that, for a big enough, 

^a^{x) < Po{x) < m^'+(a;) for all a; G K,. 

From [T], it then follows that 

ml;-{x)<p'{t,x)<ml;+{x), 

which implies the uniform boundedness of p^{t,x) in e since m^'^(x) are uniformly bounded in 
e. Then there exists a compactly supported family of probability measures nt^x on R (i.e. Young 
measures; see Tartar [22]) and a subsequence (still denoted by) p^{t, x) such that, for any continuous 
function /(p), 

f{p'{t,x)) ^ i-Kt^xJik)) ase^O. (3.5) 

On the other hand, by Section [5TT1 the sequence p'^{t,x) satisfies the entropy inequality (|2.3p for 
F^{x,p) and the steady-state solutions = m^'^. In particular, we use (j3.5p and the definition 
of the sequence F^{x,p) in (|3.ip to conclude that, as e ^ 0, the compactly supported family of 
probability measures TTt^x satisfies that, for any test function J : R+, 

J {{-Kt.x] \k - m^{x)\)dtJ + (7rt,:r;sgn (fc - m^) {F{x,k) - a)) dxJ) dxdt 

+ J \po{x)-mt{x)\j{0,x)dx>0. (3.6) 

Thus, TTt.x is a measure- valued entropy solution of (jl.ip - (jl.2p with compact support for a.e. {t, x) e 
R^ in the sense of Definition 2.2. 

3.2. Existence of entropy solutions when F is discontinuous in x. We are now ready to 
state the main theorem of this section. 

Theorem 3.1. Let F{x,p) he strictly convex or concave in p for a.e. x £ R and satisfy (Hl)- 
(H3), or let F{x,p) satisfy (H1)-(H2) and (H3'). Let po{x) G L°° . Then the sequence of entropy 
solutions p^ of the Cauchy problem p.2p (in the sense of Definition 2.1) converges to the unique 
entropy solution of the Cauchy problem (|l.ip - p.2p in the sense of Definition 2.1. 

Proof. We consider the two cases separately. 

For the case (H1)-(H2) and (H3'), that is, the flux function F is monotone in p, we apply the 
compactness framework established in Section 2 to establish the convergence. For this case, the 
existence of entropy solutions has been established in In Remark 3.1, we have shown that the 
limit of the entropy solutions p'^ is determined by ci ni6cLS lire- Veil lied- entropy solution tt^^^x- 

Then, 

by Theorems 2.1-2.2, irt^x is the Dirac measure concentrated on the unique entropy solution p{t, x) 
of (|l.ip ~ (|1.2p in the sense of Definition 2.1, which implies the whole sequence converges. 

For the case (H1)-(H3), since wc have not established the existence of an entropy solution, we 
employ the compensated compactness method to establish the convergence of the entropy solutions 
of the Cauchy problem (|3.2p . which also yields the existence of a unique entropy solution of the 
Cauchy problem pTl) -pr2 l) . 

From Remark 3.1, we have known that p'^ is uniformly bounded in L°° which implies that there 
exists a subsequence p^ converging weakly to a compactly supported family of probability measures 
vt^x on R+ such that, for any function f{p,t,x) that is continuous in p for a.e. {t,x), 

f{p^t,x),t,x) ^ {vt^x, f{k,t,x)) ase->0. (3.7) 



In particular. 



p^it^x) ^ {vt,x,k) =: p{t,x) e L' 



(3.8) 
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Our goal is to prove the strong convergence of p^{t,x) to p{t,x) a.e., equivalently, ly^x = Sp(t,x), 
which imphes that p{t, x) is an entropy solution of l|l.ip - (|1.2p . that is, p{t, x) satisfies the entropy 
inequality in Definition 12.11 

By Section 3.1, we have known that the sequence p^ exists and satisfies 

dt \ps{t, x) - p%s, y,x)\+ a^sgn {p%t, x) - p^s, y, x)) {F,{x, p'{t, x)) - 7(5, y)) ) < 

in the sense of distributions, where 

ff{s,y,x) := ■m^f^^ y^{x)sgn^{p{s,y) - pra{y)) + m:^^^^ y.^{x)sgn^{p{s,y) - Pmiy))- 

Notice that 'j{s,y) := F(y, p(s,y)) is independent of e. Thus, for fixed (s,y), we have the strong 
convergence of w^^'^ (x) to a steady-state solution rri^^^ (x) of (|l.ip - (|1.2p as e ^ 0. In particular, 

||/5^||l=o <M, M independent of e; 

and, for a.e. {s,y,x) £ x K,, 

ff{s,y,x) -> p{s,y,x) := m+ ^^^^(a;)sgn+(p(s, y) - + m:^^^^^^(a;)sgn_(p(s, y) - 

as e — *■ 0. By Schwartz's lemma, E'^ is a sequence of measures; by Murat's lemma [2(>, E'^ is 
uniformly bounded measure sequence in the measure space. This implies that 

E"" is compact in VK,;^^'^(l[l+) for any p e (1, 2). (3.9) 

On the other hand, since the vector field sequence 

i\p'it,x) - sgn{p'it,x) - m^lly^ix)) (F.ix, p%t, x)) - 7(5,2/))) 

is uniformly bounded in e for any fixed (s,y), it follows that 

E^ is bounded in VK;~^'°°(]R^). (3.10) 

With (|3.9p - p.l0p . we obtain by a compact interpolation theorem in [31[TT] that 

E'' is compact in iJ;;](]R,^). (3.11) 

On the other hand, 

dtp^ + dxFf;{x, p^) = which is automatically compact in i7;^^(R^). (3.12) 
Moreover, since /5'^(s, y, x) strongly converges a.e., then we find that, as e ^ 0, 

T]l{p'',t,x,s,y) := \p''{t,x) - p''{s,y,x)\ 

^ {i't.x{k);\k- p{s,y,x)\) 

=■■ {iyt,x;Viik,t,x,s,y)), 

9l(p^i,a;,s,y) := sgn{p''{t,x) - p''{s,y,x)) {F^{x,p^) - j{s,y)) 

^ {n,x (k); sgn (fc - p{s, y, x)) {F{x, k) - 7(5, y))) 

=■■ {i't.x;qi{k,t,x,s,y)), 

^ {'^t,xik);k) = p{t,x) 

=■■ {vt,x;m{k)), 
ql{p'{t,x),x) := F,{x,p') 

- {vt,x{k):F{x,k)) 

■= {vt.x]q2{k,x)), 



(3.13) 



and 



771 {p^{t, x), s, y, x) qi {p^{t, x), s, y, x] 
V2ip^{t,x)) q2ip%t,x),x) 



vt.. 



m{k,s,y,x) qi{k,s,y,x) 
mik) q2{k,x) 



(3.14) 
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where 

{T]i{k,t,x,s,y), qi{k,t,x,s,y)) = {\k - p{s,y,x)\, sgn {k ~ p{s , y , x))) (F{x,k) - j{s,y))), 
imik),q2{k,x)) = (fc, F{x,k)). 

Together ((XTT|) - (IXT^ with ^J^-^J^, we apply the Div-Curl lerama (see Tartar [22 and Murat 
[TO] ) to obtain 

{i^t,x; Vi{k, s, y, x)) {vt,x; qi{k, s, y, x)) 
{vt,x;ri2{k)) (vt,x;q2{k,x)) 

for all (s, y), {t, x) e R\A^ with A4 a set of measure zero in R^. Thus, we have 
{i^ty, \k - p{s,y,x)\F{x,k) - k sgn {k - p{s,y,x)) {F{x,k) - 'y{s,y))) 

= {vt.x] \k - p{s,y,x)\) {vt,x;F{x,k)) - {iyt.x,k) {i^t,x;sgn{k - p{s,y,x))) {F{x,k) - ■y{s,y))) . 
Equivalently, we have 

{i^ty, \k - p{s, y, x)\ {F{x, k) - {vt,x] F{x, k)))) 

- {vt.x\ {k - pit, x))sgn (k - p{s, y, x)) {F{x, k) ~ F{y, p{s, y)))) = 0. 

Since this is true for all (s,y) and {t,x) except on a set M of measure zero, we then choose 
(s, y) — (t, x) for {t, x) e R\M to obtain 

{'^t,x; \k ~ p{t, x)\ {F{x, k) - {yt^x\ F{x, fc)))) 

- {vt.x-, [k - p{t, x))sgn [k - p{t, x)) {F{x, k) - F{x, p{t, x)))) = 0, 

that is, 

{i^ty, \k - p{t, x)\) {Fix, p{t, x)) - {vt,x; F{x, k))) = 0. (3.15) 
There are two possibilities: 

When {vty, \k — p{t,x)\) — 0, then we have iyt,x{k) = Sp(^t.x){k)- 
When {i't,x', F{x, k)) — F{x, p{t, x)) — 0, we note that 

{iyt,x;F{x,k)) - F{x,p{t,x)) = {iyt,x;F{x,k)- F{x,p{t,x))) 

= {i^t,x;Fp{x,p){k - p) + ^Fppix,^){k ~ pf) 

= Fp{x, P) {vt,x] k- p) + ]^{vt^x\Fpp{x, £,){k - pf) 

= -{i^t,x;J eFpp{x,6p + {i-6)k)de{k^ pf). 

Since F{x, p) is strictly convex or concave in p, we conclude 

vtAk) = 5p(t.x){k) for (t, x) a.e. (3.16) 

Therefore, we have 

p'^(t,x) —> p{t,x) a.e. as £ ^ 0. 

Since the limit is unique via the uniqueness result in |T], the whole sequence p^{t,x) strongly 
converges to p{t, x) a.e. It is easy to check that p{t, x) is the unique entropy solution of the Cauchy 
problem (|l.ip - (|1.2p in the sense of Definition 2.1. □ 

Remark 3.2. The conditions on the flux function F(x,p) in Theorem 3.1 for the non-monotone 
case can be relaxed as follows: F{x,p) satisfies (H1)-(H3) and is convex or concave with 

C^{p : Fpp{x, p) = 0} = for a.e. x G R, 

where is the one-dimensional Lebesgue measure. 



m{k,s,y,x) 
V2{k) 



qi(k,s,y,x) 
q2{k,x) 
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4. Hydrodynamic Limit of a Zero Range Process with Discontinuous 

Speed-Parameter 

In Section 2, we have established a compactness framework for approximate solutions via the 
reduction of measure- valued entropy solutions of (|l.ip - (|1.2p in the sense of Definition 2.1. In this 
section we focus on a microscopic particle system for a Zero Range Process (ZRP) with discontinu- 
ous speed-parameter A(x). We apply the compactness framework to show the hydrodynamic limit 
for the particle system, when the distance between particles tend to zero, to the unique entropy 
solution of the Cauchy problem 

dtp + dAKx)Hp)) = G (4.1) 

and with initial data: 

p|t=o = Po{x) > 0, (4.2) 

where h{p) is a monotone function of p, and X{x) is continuous in a; G It with < Ai < A(x) < A2 < 
00 for some constants Ai and A2, except on a closed set J\f of measure zero. Then m+ = := TOq 
for a G [0, 00). 

Rezakhanlou in [21] first established the hydrodynamic limit of the processus des misanthropes 
(PdM) with constant speed-parameter. Covert-Rezakhanlou [TD] provided a proof of the hydro- 
dynamic limit of a PdM with nonconstant continuous speed-parameter A. In both proofs, the 
most important step is to show an entropy inequality at microscopic level, which then implies the 
(macroscopic) Kruzkov entropy inequality, when the distance between particles tends to zero, and 
thereby implies the uniqueness of limit points. In this section, we generalize this to the case when 
the speed-parameter A has jumps for the attractive Zero Range Process (ZRP). In §4.1, we analyze 
some properties of the ZRP. In §4.2, we prove the one-dimensional microscopic entropy inequality 
letting e = e(N) = N~°',a G (0, 1), for a ZRP with discontinuous speed-parameter as ^ 00. 
In §4.3, we show the existence of measure- valued solutions via the microscopic entropy inequality 
and how inequality (|2.3[) follows. 

4.1. Some properties of the microscopic interacting particle system. We consider a system 
of particles with conserved total mass and evolving on a one-dimensional lattice Z according to a 
Markovian law. With the Euler scaling factor N, the microscopic particle density is expected to 
converge to a deterministic limit as iV — > 00, which is characterized by a solution of a conservation 
law. Under the Euler scaling, represents the distance between sites. Obviously we have two 
space scales: The discrete lattice Z as embedded in R with "vertices" and m G Z. In this way, 
the distances between particles tend to zero if N increases to infinity. Sites of the microscopic scale 
Z are denoted by the letters m, v and correspond to the points , in the macroscopic scale R. 
Points of the macroscopic space scale R are denoted by the letters a;, y and correspond to the sites 
[xN], [yN] in the microscopic space scale, where [z] is the integer part of z. We denote by r]t{u) 
the number of particles at time t > at site u. Then the vector rjt = {rit{u) : u G Z) is called a 

Z 

configuration at time t with configuration space W . 

In general, the ZRP can be described as follows: Infinitely many indistinguishable particles are 
distributed on a 1-dimensional lattice. Any site of the lattice may be occupied by a finite number 
of particles. Associated to a given site u there is an exponential clock with rate Ae(-^)g (?](?/)) 
depending on the macroscopic spatial coordinates. Each time the clock rings on the site u, one 
of the particles jumps to the site v chosen with probability p{u,v). The elementary transition 
probabilities p: Z 1-^ [0, 1] are supposed to be 

(i) translation invariant: p{x, y) = p(0, y — x) —: p{y — x); 

(ii) normalized: ^yPix,y) — 1, p{x,x) — 0; 

(iii) assumed to be of finite range: p{x, y) = for \y ~ x\ sufficiently large; 

(iv) irreducible: p{0, 1) > 0. 
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Without loss of generality, we assume that YlzPi^)^ — 7 = otherwise, for 7 ^ 1, we replace 
the function h{p) by h{p)/^ in the following argument. The rate : IN ^ 11+ is a positive, 
nondecreasing function with g{Q) = 0, g{+oo) = +00, and 

4^^0 as fc->c5o. (4.3) 

With this description, the Markov process r}t is generated by 

ATLf /(^) =NY,K{^)9{'n{u))p{v - u){f{'n-n - m). (4.4) 

Here N comes from the Euler scaling factor speeding the generator, thus rjt denotes a configuration 
on which this speeded generator NL^ has acted for time t, and r/"'" represents the configuration 
7] where one particle jumped from u to v: 

{rj{ti!) li w ^ u, V, 

r]{u) - 1 if w = M, 
r]{v) + 1 a w = V. 

For any s = s{N) > and for any constant a > 0, we define a product measure given by 

'^^''^ ^■=^WJhmim^^^• -n'>.^(^M)> (4-5) 

where Z is a partition function equal to 

00 „ 

Ae(iv) ;^o(A^(#)) 9{n)\ 
Then the expected value of the occupation variable rj{u) is equal to 

Now let h be the inverse function of R to obtain 

''(^(a;^^)) = ^ Ae(;^)/i(i^£>«[»?(M)]) =a 4^ £;piv[r/(M)] =m„(-^), 

where is a steady-state solution to 

dtp + d^{X,{x)h{p))=Q. (4.7) 

Furthermore, it follows that 

Ev^iginiu))] = /i(ma(^)) . 

From now on, we set 

mK^) ^^^m^i^Mu)) := n^^e(t)M-c.(-))('?(«))- (4-8) 

u u 

The important attribute of the ZRP with nonconstant speed-parameter is that the product 
measure p^ {rj) is invariant under the generator NL^ , i.e.. 



/ 



L^{f{v))dpZAv) = 0. (4.9) 
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As a reasonable initial distribution, we choose the product measure fj,Q (rj) associated to a bounded 
density profile defined as follows: For a bounded density profile po ^ 0, the probability that 
particles at time t ^ are distributed with configuration rj is equal to 

where pu,N > is a sequence satisfying limAr^oo / \pinx],n ~ Poix)\dx — for [Nx] as the integer 
part of Nx. With this definition, we say that a sequence of probability measures fi^ is associated 
to a density profile p > if 

lim (^^(77); Itt / >^(t7)'7(w) ^ / J{x)p{x)dx\) — for every test function J. 
i^oo 'N ^ — ' A* J 



Furthermore, let 



/if (4.11) 

where — e*"^^^- is the semigroup corresponding to the generator NL^ . Then the attractiveness 
for two initial measures p^ and p^ with profiles pt and wj, respectively, implies that 



p^ <p^ ^ p"" <p^ 

is satisfied by the assumption that 17 is a nondecreasing function. Moreover, it is easy to prove 
that Ppg < Pi^g if pq < luq. It then follows by attractiveness that, for any constant a such that 
maix) > po{x), we obtain that the inequality p^ < p^^ implies 

5fMo"^<5f/xl=Ml. (4.12) 

Since our initial distribution has a bounded density profile, then the density profile remains bounded 
at later time t. 

The goal in proving the hydrodynamic limit of a ZRP is that, if we start from a configuration 
770 distributed with an initial measure Pq associated to the bounded density profile po, then the 
configuration rjt at later time t is distributed with the measure p^ defined by (j4.1ip and having 
density profile p{t, •), where p is the solution of the Cauchy problem (|4.ip - (|4.2p in the sense of 
Definition 2.1. In other words, our main theorem in this section is the following. 

Theorem 4.1 (Hydrodynamic limit of an attractive ZRP with discontinuous speed-parameter). 
Let rjt be an attractive ZRP with (|4.3p initially distributed by the measure Pq associated to a 
bounded density profile pQ : >—>■ as defined in (j4.10p . Let e ~ e(A) = N^"^ , a G (0, 1). Then, 
at later time t, 

Jim^(Mf(r7); |^E'^(^)'^*(")- / J{^)p{t,^)dx\)^Q (4.13) 

for any test function J : —> M, where p is the unique solution of the Cauchy problem (j4.ip - (j4.2p 
in the sense of Definition 2.1. 

To achieve this, we have to establish an entropy inequality in the sense of Definition 12.11 at 
microscopic level. This will be done in §4.2 by using the scaling relation e — e{N) — N~'^ , a S (0, 1). 
Associated to each configuration rjt, we may define the empirical measure viewed as a random 
measure on R by 

X^'{x)--^j^Y.^t{u)S^{x)- (4.14) 

U 

Then (xf (•), J{-)) =TrEu JifHiu), and we can rewrite by 

hm (pfir,); \{x? {■), J{-)) - / Jix)p{t, x)dx\) = 0. (4.15) 

N—nx / 
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4.2. The entropy inequality at microscopic level. The following proposition is essential to- 
wards the hydrodynamic limit. 

Proposition 4.1 (Entropy inequality at microscopic level for e = N^"' with a E (0, 1) as TV — > oo). 
Let m„ be the steady-state solutions of p.2p as defined in (|1.3p with F^[x,p) — X^{x)h(p). Let rjt 
be the ZRP generated by NL^ defined by (j4.4p and initially distributed by the measure Hq defined 
by ()4.10|) . Letrfiu) be the average density of particles in large microscopic boxes of size 2Z + 1 and 
centered at u: 

\u—v\-<l 

Then, for every test function J : — ^ 

^hm ^hrn^Mf { / ^ E (^^•^(^' J^Min) - ™^(^)| + 9. J(s, ^)| Ae(^)M^/iH) " a\)ds 
+ ^E^(0'^)|'?oW-<(^)| > -5} = 1. (4.16) 

u 

Inequality (|4.16p is the entropy inequality (|2.3p with p replaced by the average density of particles 
in the microscopic boxes of length 2Z + 1 . To prove the microscopic entropy inequality we consider 

the coupled process [rjt, it) generated by NL^, where L^ is defined by 

u,v 

+ J2piv- u)\,i^){girjiu)) .g(e(«))}+ - fiv, 0) 

u,v 

+ Y.piv- u)X,{^){giiiu)) giviu))U {f{v. C'l - fiv, 0) ■ (4.17) 

u,v 

Furthermore, denote the initial distribution of {iit,£,t) by p,Q ~ fiQ x //^e , where fiQ is the initial 
measure with density profile po defined by (I4.10p and /i^E denotes the invariant measure as defined 
in gH). 

Then, to prove Proposition [5711 it suffices to prove the following proposition. 

Proposition 4.2. Let {'r]t,£,t) be the coupled process, starting from p.Q , generated by NL^ as 
defined by (|4.17p . Let p,^ — * p.Q , where is the semigroup corresponding to the generator 
NL^ . Then, for every test function J : R+ and every e — N^'^ with a £ (0, 1), 

+ ^ E -^(0' ^) l^o(") - ^oin) I > -4 = 1. 

Recall that a microscopic entropy inequality leading to the Kruzkov entropy inequality has 
been proved in pU] for the process of PdM with nonconstant but continuous speed-parameter 
Ag. Since there does not exist an invariant product measure for a PdM in general such that 
^Ai^e ~ "^a(:^)' to replace the process by the process m'^{^), one has to apply the 

relative entropy method of Yau [23] . 



16 



GUI-QIANG CHEN, NADINE EVEN, AND CHRISTIAN KLINGENBERG 



In our case of a space-dependent ZRP, the invariant product measure is available so that we 
can approximate the steady-state solution by a process ^ distributed by the invariant measure 
/i^E for any a G (0, cx)). Then, Proposition 14. II indeed directly follows from Proposition 4.2. 

4.3. Proof of Proposition [4T2l We split the proof in three steps. 

Step 1: Lower hound for the martingale. For a test function J with compact support in M^, 
define by the martingale vanishing at time t = 0: 

u u 

Since J has compact support, then, for t large enough, 

= E -^(0' ^) - ^"Hl - fids + NL^){^ E '^(^' ^) - )'^^- 

We now calculate 

L^\n{u)-au)\ 

+{9{v{v)) ~ .g(e(«))}+( h^'^'^C^.) - (.{u)\ Hu) au)\ ) 

= E (1 - 0) {-P{v- u)K{^) \gHu)) - g{i{u))\ + p{u - t;)A,(^) |g(^(«)) - g{i{v))\ ) 

- EG«."(^'^)(p(« - ")^^(^) Is!'?!")) - - ^)Ae(^) - 9{i{v))\), (4.18) 

where Gii^t, is the indicator function that equals to 1 if 77 and ^ are not ordered, i.e., 

GuAv-X)=t{ii{u)<e,{uYj^{v)>^{v)} + t{T^{u)>^{u)-r^{v)<^{v)}. 

Notice that the second sum is nonpositive. Therefore, plugging in the last expression in the 
martingale M^ and then interchange u and v in the last term, we can bound the martingale below 

by 

u ^ u 

+ /* E (•^(^' ^) - Ji'' J^M^ - ^)(l - GuAVs,^) Ae(^) \giVs(.u)) - g{Uu))\ ds. 
Since the transition probabihty p is of finite range, i.e. p{z) = if |2;| > r for some r, then 
^) - J{s, ^)) piv ^ u) = - u)piv - u)d,Jis, ^) + Oi^). 
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With V = u + y, it then foUows that the martingale is bounded below by 

+ d^Jis, ^)A,(^)r„(^yp(y)(l - Go^y))\giVsm ~ g(6(0))|}ds 

V 

U 

Step 2: We show 

\im^E-^.[{M^f]=Q. (4.19) 

Recall that 

is a martingale vanishing at time t ^ 0, where A"^ is defined by 

u 

Then, by definition, Ep^N [iV/] = for all < s < Thus, it suffices to show that the expectation 
of the integral term of Nf' converges to zero as ~> oo. In order to prove this, we first find that, 
by careful calculation, 

NL^{A'{s, rj, Of - 2NA'{s, rj, OL^{A''{s, rj, 0) 

= J2p{w - \9iVs{v)) - giUv))\ ^(1 - G,,UVs,^s)){J{s, ^) - J{s, 



Since J is a smooth function, the first term of this expression is less 0(2^^^) for some constant 
C depending on the total initial mass and therefore converges to zero as iV ^ oo by (14. 3p . For the 
second term, we know that (J(s, f ) + J(s, f ))2 < 4 |1 J|l^, which implies 

V,W 

Then, to conclude the proof of (|4.19p . it suffices to show 

{Y.GvAVs,is)p{w - v)X,i^)\g{isiv)) - g{vs{v))\) ds] = (4.20) 

v.w 

For this, we use the martingale Mf vanishing at with J = 1, that is, 

u u ^ u 

By (|4.18p . the integral term of the martingale is equal to 

* 2 u 
^ N E NGu.vins,Mv - ")Ae(^) \9 iVs{u)) - g ds, 
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by interchanging u and v in some terms. Then we find 



'2^Gu,v{Vs:is)p{.v - u)Xe{^) \g{r]s{u)) - g {£,s{u))\ ds 



N 



Since we assumed that both marginals of /zf are bounded, (|4.20|) foUows, which leads to (|4.19p . 
With the result of Step 1 and (|4l9l) and using the Chebichev inequality, we obtain 



u ^ u 

+a.J(s,^)A,(^)r4^yp(y)(l-Go,,)(r,,e))|ff(77.(0))-5(es(0))|}ds + 0(l)<-<5} 

y 

< /If {M/ > 6} < /If {|M/| > <5} < ^i?-„[(M/)^], 
which converges to as iV ^ cxi, for all (5 > 0. 



(4.21) 



Step 3. We next use the following summation by parts formula: For any bounded function a of 
ri[-) with a(0) = and for any smooth test function J : R — + H, we obtain that, for any L > 0, 



i\\j\ 



Lip 



\u\<LN 



|u|<LAr 



\u-v\<l 



N 



). (4.22) 



Since we restrict e = 7V-'^,cr G (0, 1), then \\\e\\Lip < C/e = CN'' and ©(^IMii^) = 0{j^) 
as TV — > oo so that we can use this summation by parts formula (|4.22p to replace inequality (|4.2ip 

by 



.^-.1- E i^o(.)~eo(. 

u \z-u\<l 



\z—u\<l 

u , 1 



7V'2/ + 1 



E r4Y.yp{y)il-Go,y){Vs,Cs))\9iVsm-giCsm\ds<-5}=0. (4.23) 

2-«|<i y 

Notice that, in (|4.23p . since J is a positive function, by the triangle inequality, we can remove the 
sum inside the absolute value in the first line. Following the same argument as in [TUIEI] (also [S]), 
since we first set e — independent of Aj(a;), we can obtain the following one block estimates: 

^hm Jmi^i?^«{ ri^l-i— \Vs{z)-U^)\-\rjiiu)-es{u)\\ds}=0, (4.24) 

"'O « \u-z\<l 
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and 



lim ^lini^ E^. { r ^ E ^"1 ^ E E yP(y^^^ - Go,y){Vs,^s)) \9iVsm - giUO))\ 

-\hUiO))-hiesm\\ds}^0. (4.25) 



Combining (|4.23p with (I4.24p - (|4.25p . we complete the proof of Proposition 4.2. 

4.4. Existence of measure-valued entropy solutions. In this section, we prove that Theorem 
14.11 implies the existence of a measure- valued entropy solution associated to the configuration % ■ 
We recall the empirical measure (x) associated to a configuration rjt in (|4.14p . We define the 
Young measures associated to ry^ as follows: 

u 

which implies {tt^'^;J) ~ ;^ X^u "^(w' '^t(^)) ^^^^ ^ Co{R x R+). If E is the configuration 
space, then these two measures are finite positive measures on E and, for any J G Co(]R), they are 
related by the formula 

(7rf-';fcJ(x))«(xf(-);J(-)>- (4.27) 

Notice that, since there are jumps, the probability measure fi^ defined by (|4.1ip must be defined 
on the Skorohod space D[(0, oo), E], which is the space of right continuous functions with left limits 
taking values in E. Then, using the one to one correspondence between the configuration rjt and 
the empirical measure Xt'(')i the law of with respect to will give us a probability measure 
on the Skorohod space D[{0, oo), 7U+(1R)], for the space M+(R) of finite positive measures on 
R endowed with the weak topology. 

In the same way, we can associate a probability measure Q^'' on the space D[{0, oo), A1+(R^)]. 
With these definitions, we can state the main theorem of this section as follows. 

Theorem 4.2 (Law of large numbers for the Young measures). Let (/i^)7v>i be a sequence of 
probability measures, as defined by (I4.10[) . associated to a bounded density profile po ■ R ^ R+- 
Then the sequence Q^'' converges, as N ^ oo first and I — > oo second, to the probability measure 
Q concentrated on the measure-valued entropy solution 'Kt,x in the sense of Definition 2.2. 

Proof. In order to be allowed to take the limit points Q and Q of and Q^'' respectively, we 
must know that the sequences are tight (weakly relatively compact). If Q^'' is weakly relatively 
compact, we can take as a limit point if — > oo for each I. Denote by Q a limit point of 
Q^'i \{ N ^ oo first and / oo second. Therefore, the proof consists in two main steps: The 
first is to show that Q^'^ is weakly relatively compact and the second is to show the uniqueness 
of limit points. The key point in the proof is that these can be achieved independent of the 
choice of mollification of the discontinuous speed-parameter A with our choice of the notion of 
measure- valued entropy solutions. 

These can be achieved by following exactly the standard argument in [10l[2Tl[T6] since it requires 
only the uniform boundedness of A^ in the proof. That is, we can conclude the following: Let 
be a measure defined by (|4.1ip . Then 

(i) The sequence defined above is tight in Z?[(0, oo), A1 + (R)] and all its limit points Q 
are concentrated on weakly continuous paths x(t, •); 

(ii) Similarly, the sequence Q^'' is tight in i5[(0, oo), A^_|_(R x R+)] and all its limit points Q 
are concentrated on weakly continuous paths 7r(t, •, •); 
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(iii) For every t>0, 7r(i, a;, k) := TTt{x, k) is absolutely continuous with respect to the Lebesgue 
measure on R, Q a.s.. That is, Q a.s. 

T:t{x,k) =T:t,x{k)®dx; (4.28) 

(iv) For every t € [0, T], 'Kt,x{k) is compactly supported, that is, there exists fcp > such that 

^M([o,fco]") = V(i,x) e [o,r] X R 

(v) 'Kt,x is a measure- valued entropy solution in the sense of Definition 2.2 for any a G [Afg, oo), 
i.e., 

dt {TTty, \k - ma{x)\) + a, {TTt,x; \h{k)X{x) - a|) < (4.29) 

on in the sense of distributions for any a £ [Mq, oo) or a e (— oo, Mq]. 

The last result follows from the entropy inequality at microscopic level in Theorem 14. II Indeed, 
in terms of the Young measures, the expression (|4.16p of Proposition 4.1: 

+d,H{t, ^)\Xi^)h{4{u)) - a\ }dt >-5} = l 



can be restated as 



\N,l \ 



lim lim Q'^''\ / U'Kt{x,k);\k-ma{x)\dtH(t,x)) 

l—^OO N—>-oo 



+ {TTtix,k);\X{x)h{k) ~ a\d^H{t,x)))dt > -S^ = 1, 

for every smooth function H : (0, T) x R — )■ R_|_ with compact support, any a G [Mq, oo) or 
a e (— oo, Mq], and any S > 0. Since Q is a weak limit point concentrated on absolutely continuous 
measures and since we already proved that nt.x is concentrated on a compact set (and therefore 
the integrand is a bounded function), we obtain from the last expression that 

q{ f n {nt.x; \k - m„(x)|) dtH{t,x) + (tt*,,; \X{x)h{k) - a\) d^H{t, a;)) dxdt > -6^ = 1. 



Letting S ^ 0, we have that Q a.s. (|4.29p holds on (0, T) x R in the sense of distributions for 
every a £ [0, oo). This proves the uniqueness of Q and thereby concludes the proof of Proposition 
im □ 



Then Theorem 14.11 follows immediately from this result since the measure- valued entropy solu- 
tion reduces to the Dirac mass concentrated on the unique entropy solution p{t, x) of (|4.ip - (|4.2p 
as we noticed in §3.2. 
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